Surface-emitting  circular  DFB,  disk-,  and  ring- 
Bragg  resonator  lasers  with  chirped  gratings. 
Ill:  gain  saturation  effects  and 
above-threshold  analysis 

Xiankai  Sun*  and  Amnon  Yariv 

Department  of  Applied  Physics,  MC  128-95,  California  Institute  of  Technology,  Pasadena,  CA  91125,  USA 
Corresponding  author:  xksun  @  caltech.  edu 


Abstract:  As  Part  III  of  this  series,  this  paper  focuses  on  an  above¬ 
threshold  modal  analysis  which  includes  gain  saturation  effects  in  the 
surface-emitting  chirped  circular  grating  lasers.  We  derive  an  exact  energy 
relation  which  states  that,  in  steady  state,  the  net  power  generated  in  the 
gain  medium  is  equal  to  the  sum  of  peripheral  leakage  power  and  vertical 
emission  power.  This  relation  is  particularly  useful  in  checking  the  accuracy 
of  numerical  mode  solving.  Numerical  simulations  demonstrate  the 
dependence  of  required  pump  level  on  the  vertical  emission  power  and 
compare  the  laser  threshold  and  energy  conversion  efficiency  under 
uniform,  Gaussian,  and  annular  pump  profiles.  A  larger  overlap  between  the 
pump  profile  and  modal  intensity  distribution  leads  to  a  lower  threshold  and 
a  higher  energy  conversion  efficiency.  Finally  the  dependence  of  required 
pump  level  on  device  sizes  offers  us  new  design  guidelines  of  these  lasers 
for  single-mode,  high-efficiency,  high-power  applications. 
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waves;  (230.1480)  Bragg  reflectors;  (050.2770)  Gratings;  (140.5960)  Semiconductor  lasers; 
(130.0130)  Integrated  optics. 
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1.  Introduction 

Surface-emitting  circular  grating  lasers  are  considered  a  promising  next-generation  on-chip 
light  source  for  optical  communications  which  may  substitute  the  prevailing  VCSELs  due  to 
their  circularly-symmetric,  large-emission-aperture,  narrow-divergence  laser  emission.  There 
have  been  intensive  research  activities  since  early  1990s.  Erdogan  and  Hall  were  the  first  to 
analyze  their  modal  behavior  with  a  scalar  coupled-mode  theory  [1],  Wu  et  al.  were  the  first 
to  experimentally  realize  such  lasers  in  semiconductors  [2].  More  recently,  organic  polymers 
are  also  used  as  the  gain  medium  for  these  lasers  due  to  their  low  fabrication  cost  [3-5]. 

The  circular  gratings  in  most  of  the  work  are  designed  radially  periodic.  In  2003  we 
proposed  using  radially  chirped  gratings  to  achieve  optimal  interaction  with  the  optical  fields, 
since  the  eigenmodes  of  the  wave  equation  in  cylindrical  coordinates  are  Hankel  functions  [6]. 
As  the  grating  designs  follow  the  phases  of  Hankel  functions,  such  circular  grating  lasers 
usually  take  three  configurations  shown  in  Fig.  1(a)  circular  DFB  laser,  (b)  disk  Bragg 
resonator  laser,  and  (c)  ring  Bragg  resonator  laser.  Including  a  second-order  Fourier 
component,  the  gratings  are  able  to  provide  in-plane  feedback  as  well  as  couple  laser  emission 
out  of  the  resonator  plane. 


Fig.  1.  Surface-emitting  chirped  circular  grating  lasers:  (a)  circular  DFB  laser;  (b)  disk  Bragg 
resonator  laser;  (c)  ring  Bragg  resonator  laser.  Laser  emission  is  coupled  out  of  the  resonator 
plane  in  vertical  direction  via  the  gratings. 


In  previous  parts  of  this  paper  series  [7,  8],  we  have  solved  the  modes  and  compared  the 
modal  properties,  by  both  analytical  and  numerical  approaches,  of  the  three  types  of  lasers  at 
threshold.  Gain  saturation  effects  in  above-threshold  operation  were  previously  studied  in 
both  linear  and  circular  DFB  and  DBR  lasers  [4,  9-12].  This  paper  presents  an  above¬ 
threshold  analysis  for  the  three  types  of  radially  chirped  circular  grating  lasers  to  investigate 
their  energy  conversion  efficiency,  intensity-dependent  modal  selectivity,  and  optimal  design. 

2.  Energy  relation 

The  basic  equations  that  govern  the  modal  behaviors  in  such  Hankel-phased  active  circular 
grating  structures  have  been  derived  in  [13]: 

=  u(x')  ■  A(x')  -  v  ■  B(x)  ■  e2lSx ,  (la) 

dr 

=  —u  (x)  •  B(x)  +  v  ■  A(x)  ■  e~2lSx ,  (lb) 

dr 

where 

A(.r)  and  B(x):  amplitudes  of  the  in-plane  outward  and  inward  propagating  cylindrical 
waves; 

r  =  Pp\  normalized  radial  coordinate  with  fi  being  the  in-plane  propagation  constant; 

3  =  (/ff design ~P)IP'-  frequency  detuning  factor,  representing  the  relative  frequency  shift  from 
the  optimal  coupling  design; 

u(x)  =  -  hi  and  =  hi  +  ih2\ 

hi  =  liir  +  ih],:  radiation  coupling  coefficient,  representing  the  effects  of  vertical  emission; 
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h2 :  feedback  coupling  coefficient  (can  always  be  chosen  real); 
g,\(x)  =  g(x)  -  a:  space-dependent  net  gain  coefficient; 

a:  nonsaturable  internal  loss,  including  absorption  and  nonradiative  scattering  losses; 
g(x)  =  go(x)/[  1  +  /(x)//saJ:  intensity-dependent  saturated  gain  profile; 
go(x):  unsaturated  gain  profile;  and 

/(x)//sat:  field  intensity  distribution  in  units  of  saturation  intensity. 

All  the  three  types  of  circular  grating  lasers  share  the  same  boundary  conditions  [7]:  (i)  A(0)  = 
B( 0);  (ii)  B(xh)  =  0;  (iii)  A (x)  and  B{x)  continuous  for  0  <  x  <  xh.  Once  Eqs.  (la)  and  (lb)  are 
solved  with  these  boundary  conditions,  the  in-plane  electric  field  E(x)  =  A(x)Hm  +  B(x)H<2> . 

Assuming  a  linear  pump-gain  relationship  above  transparency,  the  unsaturated  gain  g0(x) 
follows  the  profile  of  pump  intensity  /pump(x),  and  we  may  define  an  experimentally 
measurable  quantity:  pump  level  Ppump  =  I  /pumpW  '  2 up  ■  dp  =  P0  \  ga(x')  ■  x  ■  dx,  where  P0 
having  a  power  unit  is  an  experimental  parameter  determined  by  pumping  scheme  (electrical 
or  optical)  and  specific  setup  configuration.  For  some  simple  g0  profiles,  Ppump  has  analytical 
expressions  [8],  otherwise,  numerical  integration  always  remains  a  resort. 

Similar  to  the  procedure  in  [9],  multiplying  Eq.  (la)  by  A  and  Eq.  (lb)  by  B  ,  then  adding 
each  equation  to  its  complex  conjugate,  one  obtains 

-  =  2  ( gA  -  \r )  |  A|2  -  v  •  A’ B  ■  e2,Sx  -  v*  •  AB"  ■  e~2iSx ,  (2a) 

=  -2 (gA- hlr )\B\2  +v-ABt  ■  e-2iSx  +  v*  •  A* B ■  e2iSx.  (2b) 

Subtracting  Eq.  (2b)  from  Eq.  (2a)  yields 

-^-(|a|2  -  |s|2 )  =  2gA  (|a|2  +  \b\2  )  -  2 hlr  | Ae~iSx  +  BeiSx |2 .  (3) 


dr 

d  Lb 


Integrating  Eq.  (3)  from  x  =  0  to  x  =  xb  and  applying  A(0)  =  B ( 0 )  and  B(xh)  =  0  lead  to 


|A(xfc)|2  + 

peripheral  leakage 


2 hlr  f  Ae~iSx  +  BeiSx  2  dx  =  2P  gA  (|a|2  +  |S|2)dr, 

J  grating  J  0  V  / 

v _ _ _ /  v _ w _ / 

v  v 

vertical  laser  emission  power  generated  in  the  gain  medium 


(4) 


which  is  the  energy  conservation  theorem  for  the  surface-emitting  circular  grating  lasers.  This 
equation  is  exact  and  states  that,  in  steady  state,  the  net  power  generated  in  the  gain  medium 
equals  the  sum  of  peripheral  leakage  power  and  vertical  emission  power.  Due  to  its  exactness, 
this  relation  is  particularly  useful  in  monitoring  the  accuracy  of  numerical  mode  solving. 

It  should  be  noted  that  the  amplitudes  A  and  B  are  dimensionless,  and  they  can  be  scaled 
up  or  down  arbitrarily.  A  judicious  choice  would  be  to  normalize  them  such  that  the  resulting 
intensity  is  normalized  to  the  saturation  intensity  7sat.  Since  I(x)  =  cne0\E\2/2  ( c ,  the  speed  of 
light;  n,  transverse  effective  index;  So,  the  vacuum  permittivity),  if  £sat  is  similarly  defined  as 
(2/sat/cnco)  ,  then  the  peripheral  leakage  power  in  real  dimensions  will  be 

\A(xb)H,x,(xb)\  E2dl  ■  2jixhD j p  ~  A(xh) 2  ■  E2M  AD/ p ,  where  D  is  the  thickness  of  the  laser 

resonator.  A  comparison  of  this  expression  with  the  first  term  on  the  left-hand  side  of  Eq.  (4) 
indicates  that  all  the  power  terms  are  in  units  of  the  saturation  power  p  =  E2d  AD/ (3 . 


3.  Numerical  method 

When  solving  the  modes  at  threshold  with  uniform  gain  (or  pump)  distribution  across  the 
device,  gain  saturation  effects  can  be  ignored  and  Eqs.  (la)  and  (lb)  have  analytical  solutions 
[7],  In  case  of  using  a  nonuniform  pump  profile  and/or  taking  into  account  the  gain  saturation 
effects  in  above-threshold  operation,  Eqs.  (la)  and  (lb)  have  to  be  numerically  solved: 
starting  with  an  amplitude  set  [A  B\  =  A(0)[1  1]  at  the  center,  these  equations  are  integrated 
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along  x  to  the  exterior  boundary  xh.  The  absolute  value  of  B(xh)  marks  a  contour  map  in  the 
two-dimensional  plane  of  g0  and  d.  Each  minimum  point  in  the  contour  map  represents  a 
mode  with  corresponding  g0  and  S.  Substituting  the  modal  g0,  S,  A(x),  and  B(x)  into  Eq.  (4), 
one  immediately  finds  the  mode-solving  accuracy  by  comparing  the  left-hand  and  right-hand 
sides  of  the  equation.  The  second  term  on  the  left-hand  side  of  Eq.  (4)  is  the  surface  emission 
power  Pem  in  units  of  the  saturation  power  Esat.  Inserting  the  modal  g0  into  the  expression  for 
Epump  gives  the  required  pump  level. 

During  the  numerical  integration,  an  expression  for  the  radial  intensity  distribution  I(x)  is 
required.  The  expression  for  the  electric  field  E( x)  =  AH ™  +  BH{2)  leads  straightforwardly  to 


7(x)  =  \E{xf  =  \AHa)\2  +  \BHa)\2  +  A*BHa)*Ha)  +  AB*HmH{2)* . 

v  '  |  v  ■/\  \  m  m  mm  mm 

interference  terms 


(5) 


In  a  gain  medium,  the  motional  nature  of  the  electrons  prevents  them  from  sensing  the  fast- 
oscillating  interference  terms  in  Eq.  (5).  Therefore,  with  large-radius  approximations  of  H^1 

and  H<2) ,  Kasunic  et  al.  dropped  these  interference  terms  and  simplified  the  expression  [10]. 

However,  the  remaining  first  two  terms  blow  up  at  x  =  0,  which  results  in  a  dubious 
representation  of  I(x)  near  the  center.  Here  we  suggest  approximating  7(x)  piecewise: 

*>-v 

The  criterions  for  choosing  the  joint  point  xeq  are:  first,  to  ensure  the  continuity  of  7(x),  xeq  has 
to  be  where  the  interference  terms  vanish,  i.e.,  at  a  zero  of  the  interference  terms;  second,  xeq 
has  to  precede  the  start  of  the  fast  oscillation  of  7(x);  third,  the  choice  of  xeq  has  to  validate  the 
approximation  for  7(x)  at  x  >  xeq.  Based  on  the  modal  profiles,  xeq  is  chosen  to  be  the  second 
zero  of  the  interference  terms,  which  is  1.76. 

4.  Numerical  results  and  discussions 


7(x)  = 


A77™  +  BH 


(2) 


2(|A|2  +  |B|2)//;rx, 


The  laser  structural  parameters  are  assumed  to  be  the  same  as  those  in  [7]  and  [8].  In  short,  in 
grating  regions  hx  =  0.0072  +  0.0108/  and  h2  =  0.0601  while  in  no-grating  regions  both  vanish. 
For  comparison  purposes,  all  the  three  types  of  lasers  share  an  xb.  In  addition,  for  disk  Bragg 
resonator  laser  x0  =  x,/2,  and  for  ring  Bragg  resonator  laser  xL  =  x,/2  -  n  and  xR  =  xh!2  +  n.  The 
only  parameter  related  to  material  is  the  nonsaturable  internal  loss  a,  which  is  assumed  to  be 
0.2  x  10”3  (normalized  by  fi)  for  typical  III-V  quantum  well  lasers. 

4.1.  Gain  saturation  effects  on  the  laser  modes 


(a)  Circular  DFB  laser 


(b)  Disk  Bragg  resonator  laser  (c)  Ring  Bragg  resonator  laser 


Fig.  2.  Pump  level  Ppa mp  of  the  5  lowest  order  modes,  under  uniform  pump  profile,  of  (a) 
circular  DFB,  (b)  disk-,  and  (c)  ring-  Bragg  resonator  lasers  with  surface  emission  power  Pem  = 
0.01,  1,  and  10  (in  units  of  saturation  power  PSat)- 
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To  demonstrate  the  gain  saturation  effects  on  the  laser  modes,  a  typical  device  size  pb  =  17.5 
pm  ( Xb  =  ftph  ~  200)  is  assumed  for  all  the  circular  DFB,  disk-,  and  ring-  Bragg  resonator 
lasers.  Figure  2  displays  pump  level  Ppump  of  the  5  lowest  order  modes,  under  uniform  pump 
profile,  of  these  three  types  of  lasers  with  surface  emission  power  Pem  =  0.01,  1,  and  10  (in 
units  of  saturation  power  Psat).  The  features  are  summarized  as  follows:  (i)  At  low  surface 
emission  Pem  =  0.01  and  1,  disk  Bragg  resonator  laser  always  possesses  the  lowest  Ppump  «rit  of 
the  three  types,  while  at  high  surface  emission  Pem  =10,  Ppump  of  Mode  1  of  disk  Bragg 
resonator  laser  is  greatly  enhanced  and  even  higher  than  that  of  circular  DFB  and  ring  Bragg 
resonator  lasers,  (ii)  Disk  Bragg  resonator  laser  has  smaller  discrimination  between  the  modes 
compared  with  that  of  circular  DFB  and  ring  Bragg  resonator  lasers,  and  as  a  result  the 
difference  in  gain  saturation  levels  can  easily  change  the  order  of  the  modes:  Mode  2  exhibits 
even  lower  Ppump  than  Mode  1  at  Pem  =  10.  This  mode  transition  behavior  will  be  illustrated 
more  clearly  in  the  device-size-dependent  Ppump  curves  shown  in  Fig.  3. 

4.2.  Above-threshold  behavior:  comparison  of  threshold  pump  level  and  energy  conversion 
efficiency  with  different  pump  profiles 

By  varying  the  boundary  condition  A(0)  in  the  integration  process,  we  are  able  to  get  the 
(Ppump,  Pem)  pairs  which  basically  form  the  typical  input-output  relation  for  a  laser  mode.  The 
laser  threshold  Pth  is  the  pump  level  at  the  onset  of  surface  laser  emission.  The  energy 
conversion  efficiency  (or  external  differential  efficiency)  ;/ex  is  defined  as  the  slope 
dPem/dPpump  of  the  linear  fit  of  the  simulated  data  points  up  to  Pem  =  10.  To  compare  the 
effects  of  different  pump  profiles,  Pth  and  ;/ex  of  Mode  1  of  the  three  types  of  lasers  were 
calculated  under  uniform,  Gaussian,  and  annular  pump  profiles.  The  Gaussian  profile  is 
assumed  to  follow  g0(x)  =  g0  exp^-jr2/^  J  with  WP  =  xf 2  =  100.  The  annular  profile  is 

assumed  to  follow  g0U)  =  g0exp[-U-xp)2/w^]  +  g0exp[-(x  +  xi,)7w^]  with  x„  =  xb/2  = 
100  and  wp  =  xfA  =  50.  The  results  are  listed  in  Table  1. 


Table  1.  Threshold  pump  level  Pth  (in  units  of  Po)  and  energy  conversion  efficiency  //es  (in  units  of  Psat/Po)  of 
circular  DFB,  disk-,  and  ring-  Bragg  resonator  lasers  under  different  pump  profiles 


Pump  profile 

Circular  DFB  laser 

Disk  Bragg  resonator  laser 

Ring  Bragg  resonator  laser 

Ah 

7ex 

P,h 

*7ex 

Pth 

7ex 

Uniform 

9.760 

0.7369 

6.565 

0.4374 

13.162 

0.9278 

Gaussian 

5.967 

0.9961 

2.373 

0.8741 

8.570 

1.379 

Annular 

6.382 

0.9742 

5.855 

0.7358 

7.010 

1.500 

Clearly  a  low  Pth  and  a  high  ;/ex  always  come  together.  The  lowest  /J,|l  and  the  highest  ;/ex 
are  achieved  with  the  Gaussian  pump  for  circular  DFB  and  disk  Bragg  resonator  lasers  and 
with  the  annular  pump  for  ring  Bragg  resonator  laser.  Actually  Pth  and  t/ex  are  both  related  to 
the  overlap  between  pump  profile  and  modal  intensity  distribution  [8].  A  larger  overlap  results 
in  more  efficient  energy  conversion  in  the  gain  medium  which  consequently  leads  to  a  lower 
P,h  and  a  higher  ;/ex. 

4.3.  Revisit  of  modal  selectivity  and  optimal  design 

Figure  3  shows  device-size-dependent  pump  level  Ppump  and  frequency  detuning  factor  3  of 
the  3  lowest  order  modes,  under  uniform  pump  profile,  of  the  three  types  of  lasers.  The  modes 
are  numbered  in  accord  with  those  in  Fig.  2.  For  both  Ppump  and  3,  dashed  lines  mark  their 
values  obtained  at  threshold  and  solid  lines  at  Pem  =  10. 

Examined  from  the  plots,  at  high  emission  level  Pem  =10,  circular  DFB  and  ring  Bragg 
resonator  lasers  still  possess  large  discrimination  between  the  modes,  which  ensures  them  a 
large  single-mode  range  of  at  least  50-250.  Additionally,  we  identify  low-pump  ranges  for 
their  Mode  1  at  Pem  =  10,  which  are  100-160  for  circular  DFB  laser  and  80-130  for  ring 
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Bragg  resonator  laser.  The  existence  of  a  low-pump  range  is  a  result  of  competition  between 
the  pumped  area  and  the  required  gain  level:  although  larger  devices  require  a  larger  area  to 
be  pumped,  their  longer  radial  Bragg  gratings  reduce  the  needed  gain  because  of  stronger 
reflection  of  the  optical  fields  from  the  gratings.  The  low-pump  range  is  an  important  factor  in 
designing  high-efficiency,  high-power  lasers. 

The  Ppump  for  disk  Bragg  resonator  laser  exhibits  interesting  behaviors:  (i)  at  xb  =  200,  the 
order  of  Modes  1  and  2  exchanges  from  at  threshold  to  at  high  surface  emission  level  due  to 
the  participation  of  gain  saturation  effects;  (ii)  the  single-mode  range  (for  Mode  2)  shifts  from 
60-140  at  threshold  to  90-175  at  Pem  =  10.  Therefore  the  single-mode  range  for  engineering 
should  be  the  overlap  of  these  two  sets,  which  is  90-140. 

On  the  other  hand,  all  the  frequency  detuning  factors  d  are  almost  unaffected  by  the 
surface  emission  level  (as  seen  from  the  overlap  of  dashed  and  solid  lines  in  the  lower  plots), 
since  d  is  an  intrinsic  property  of  a  laser  mode. 


(a)  Circular  DFB  laser  (b)  Disk  Bragg  resonator  laser  (c)  Ring  Bragg  resonator  laser 


Fig.  3.  Device-size-dependent  pump  level  Ppump  and  frequency  detuning  factor  S  of  the  3  lowest 
order  modes,  under  uniform  pump  profile,  of  (a)  circular  DFB,  (b)  disk-,  and  (c)  ring-  Bragg 
resonator  lasers,  xb  is  the  exterior  boundary  radius  for  all  types  of  lasers.  The  inner  disk  radius 
xo  of  disk  Bragg  resonator  laser  is  set  to  be  xJ2.  The  inner  and  outer  edges  of  the  annular  defect 
of  ring  Bragg  resonator  laser  are  set  to  be  xL  =  xJ2  —  n  and  xR  =  xJ2  +  it.  The  modes  are  labeled 
in  accord  with  those  shown  in  Fig.  2.  Dashed  lines  mark  the  values  obtained  at  threshold  and 
solid  lines  at  Pcm  =10. 


5.  Conclusions 

An  above-threshold  analysis,  including  gain  saturation  effects,  on  the  modal  properties  of 
surface-emitting  chirped-grating  circular  DFB,  disk-,  and  ring-  Bragg  resonator  lasers  has 
been  performed.  An  exact  energy  relation  in  such  surface-emitting  lasers  was  derived,  which 
is  interpreted  as  the  energy  conservation  theorem:  in  steady  state,  the  net  power  generated  in 
the  gain  medium  is  equal  to  the  sum  of  peripheral  leakage  power  and  vertical  emission  power. 

Numerical  simulations  demonstrated  the  dependence  of  required  pump  level  on  the 
vertical  emission  power  for  each  laser  mode.  The  threshold  pump  level  /Jth  and  energy 
conversion  efficiency  r/ex  were  compared  under  three  pump  profiles  -  uniform,  Gaussian,  and 
annular.  A  larger  overlap  between  the  pump  profile  and  modal  intensity  distribution  leads  to  a 
lower  Pth  and  a  higher  //ex. 
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Finally,  as  shown  in  the  evolution  of  modal  pump  levels  with  the  device  size,  disk  Bragg 
resonator  laser  exhibits  mode  transition  behavior  which  modifies  its  single-mode  range,  while 
circular  DFB  and  ring  Bragg  resonator  lasers  find  each  low-pump  range.  These  numerical 
results  provide  us  useful  information  in  designing  such  lasers  for  single-mode,  high- 
efficiency,  high-power  applications. 
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